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Abstract 

The use of human whole blood for transcrlptomic analysis has potential advantages over the use of isolated immune cells 
for studying the transcriptional response to pathogens and their products. Whole blood stimulation can be carried out in a 
laboratory without the expertise or equipment to isolate immune cells from blood, with the added advantage of being able 
to undertake experiments using very small volumes of blood. Toll like receptors (TLRs) are a family of pattern recognition 
receptors which recognise highly conserved microbial products. Using the TLR2 ligand {Pam3CSK4) and the TLR4 ligand 
(LPS), human whole blood was stimulated for 0, 1, 3, 6, 12 or 24 hours at which times mRNA was isolated and a comparative 
microarray was undertaken. A common NFkB transcriptional programme was identified following both TLR2 and TLR4 
ligation which peaked at between 3 to 6 hours including upregulation of many of the NFkB family members. In contrast an 
interferon transcriptional response was observed following TLR4 but not TLR2 ligation as early as 1 hour post stimulation 
and peaking at 6 hours. These results recapitulate the findings observed in previously published studies using isolated 
murine and human myeloid cells indicating that in vitro stimulated human whole blood can be used to interrogate the early 
transcriptional kinetic response of innate cells to TLR ligands. Our study demonstrates that a transcriptomic analysis of 
mRNA isolated from human whole blood can delineate both the temporal response and the key transcriptional differences 
following TLR2 and TLR4 ligation. 
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introduction 

Microarray analysis is increasingly being used to advance our 
understanding of the complex transcriptional responses generated 
downstream of an experimental perturbation or during a disease 
state [1,2]. Using modern microarray platforms it is possible to 
measure the expression of over 40,000 mRNA transcripts, 
encompassing all of the known functional human genome. By 
sequentially sampling over time, the temporal dynamics of the 
transcriptional response to a given event can be delineated and 
depending on the scale of the response many hundreds or 
thousands of significantly dififerentiaUy expressed genes may be 
identified. To better understand and interpret this complex data 
bioinformatics tools have been developed that take advantage of 
known biological relationships which may influence gene expres- 
sion. Using these tools it is possible to use an unbiased 
methodology to determine distinct classes of differentially regulat- 
ed genes as well as potentially important transcriptional regulatory 



pathways or networks which may change when exposed to a given 
stimulus [2^]. 

ToU-like receptors (TLR) are a family of pattern recognition 
receptors (PRR) which recognise highly conserved microbial 
products. In humans 10 functional TLRs have been described 
including TLR2 and TLR4 which are expressed on the cell surface 
and recognise conserved bacterial products [5] . TLR2 recognises 
lipoteichoic acids of Gram-positive bacteria and bacterial lipo- 
proteins whereas TLR4 recognises lipopolysaccharides (LPS) of 
Gram-negative bacteria [5,6]. TLR2, TLR3 and TLR4 are 
expressed and functional in a wide variety of cells found in human 
whole blood including dendritic cells (DC) and monocytes [6]. 
Early in the immune response to a pathogen, ligation of TLRs 
induces gene transcription leading to inflammation, tissue repair 
and the initiation of adaptive immunity [6,7]. All TLRs utilise the 
MyD88 adaptor molecule except for TLRS which uses the TRIF- 
TRAM adaptor molecules only. TLR2 and TLR4 both use the 
MyD88-TIRAP adaptor molecules. Additionally TLR4 also uses 
the TRIF-TRAM adaptor molecules [5-7]. Signalling via the 
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adaptor molecule MyD88 results in activation of transcription 
factors such as NFkB and AP-1 [5-7], whereas signalling via the 
TRIF-TRAM adaptors leads to activation of the interferon 

regulator)' transcription factors (IRF) including IRF3 [6,8]. IRF3 
is constitutively expressed and its activation by ligation of TLR4 or 
TLR3 results in to induction of IFN|3 (amongst other cytokines), 
which via the IFNaP receptor leads to positive feedback regulation 
of type 1 interferon (IFN) genes including the type 1 IFN inducible 
transcription factor IRF7 [8] . There is also a differential temporal 
transcriptional response following NFkB activation in macrophag- 
es, which can be broadly categorised into three phases: early 
primary response, late primary response and secondary response 
genes. This differential activation of sets of genes depends on their 
chromatin status and the potential need for remodeUing of the 
target genes which can lead to the differential kinetics of induction 
[7,9] . Much of the knowledge we have gained of these signalling 
mechanisms are derived from studying TLR ligation of mono- 
cytes, macrophages or DCs from mice or humans. 

A study collectively analysing data from multiple published 
papers from both human and murine macrophage transcriptional 
studies identified a group of genes upregulated following both 
TLR2 and TLR4 stimulation, which were predicted to be 
regulated by NFkB. Additionally a separate IFN-sensitive response 
element (ISRE) regulated set of genes expression was seen to be 
upregulated in the TLR4 (and TLR3) stimulations but not the 
TLR2 stimulation [10]. Importantly the temporal kinetics of these 
two groups of genes differed, with the genes predicted to be 
regulated by NFkB peaking in expression earlier than the 
predicted ISRE regulated genes. In addition, genes thought to 
be regulated by both ISRE and NFkB had a greater magnitude of 
induction though similar temporal kinetics to the genes regulated 
only by ISRE [10]. However this study had limited time points 
available for the TLR2 stimulations and so the majority of the 
temporal analysis was undertaken using LPS stimulated samples. 
Another comparative study including TLR2, TLR3 and TLR4 
stimulation of murine primary DCs identified an "inflammatory 
programme" mediated by TLR2 ligation and an "anti-viral 
programme" mediated by TLR3 Ugation. TLR4, owing to its use 
of both MyD88-TIRAP and TRIF-TRAM adaptor molecules, 
result(;d in the activation of both "inflammatory" and "anti-xiral" 
programmes [11]. The "inflammatory" programme was enriched 
for genes predicted to be regulated by NFkB (RELA, NFkBI, 
NFkBIZ) and the "anti-viral" programme was enriched for genes 
predicted to be regulated by the transcription factors STATl, 
STAT2, STAT4, IRF8 and IRF9 [11]. Study of the temporal 
transcriptomic response following TLR4 ligation in human and 
murine macrophages has also been used to specifically analyse the 
expression of transcription factors [12-14]. For example in LPS 
stimulated human macrophages the majority of transcription 
factors whose expression was seen to change had done so by 2 
hours. Groups of transcription factors peaked in expression at 
different times, coinciding with the transcriptional peaks of the 
genes they were predicted to regulate following TLR stimulation 
in both murine and human macrophages [12,15]. However studies 
in isolated cells may not reflect the overall host response to TLR 
ligation, since interaction will occur between different cell types, 
leading to a complex interplay of autocrine and paracrine 
signalling events resulting in differentiation, proliferation, cell 
trafficking and further chemokine/cytokine production and 
feedback loops for positive and negative effects on gene regulation. 

Transcriptomic studies following in vivo TLR4 ligand adminis- 
tration have also previously been undertaken [16-18]. From 
human whole blood leukocytes obtained 0, 2, 4, 6, 9 and 24 hours 
following in vivo administration of LPS, it was shown that mRNA 



expression of proinflammatory chemokines and cytokines (TNF, 
ILIA, ILIB, CXCLl, CXCL2, CCL2, CXCL8 and CXCLIO) 
peaked at 2 to 4 hours after LPS administration, whereas the 
cytokine ILIO was maximal at 6 hours. In this study the 
transcription factors NFKBI, NFKB2, RELA and RELB were 
significantly expressed and seen to peak after the cytokines and 
chemokines. The peak time for transcription factors including the 
STAT (signal transducer and activator of transcription) and IRF 
g(;nes was at 4 to 6 hours [16,19]. Analysis of mRNA isolated from 
circulating human neutrophils after 0, 2, 4 and 6 hours following 
LPS administration revealed significant upregulation of the TNF 
signalling pathway and NFkB genes such as NFkBI and NFkB2 
by 2 hours [18]. Although it is possible to undertake certain 
human in vivo experiments, studies are limited by ethical and 
practical considerations. 

Whole blood comprises cells of both the innate and adaptive 
immune system. Therefore the use of whole blood for in vitro 
studies has the potential advantage over isolated cell populations as 
these different components may havi; a diflFcrcntial response to 
stimulation. Autocrine and paracrine signalling between the 
differing cell populations may result in a response of the whole 
system that potentially better reflects the in-vivo response. 
Additionally whole blood potentially has advantages over PBMC, 
DCs or monocyte derived macrophages since it can be used 
in situations where it is not possible to obtain large volumes of 
blood to derive the isolated cell populations. Previously, In-vitro 
human whole blood has been used as a model to study TLR 
ligation, predominantiy with measurement of specific cytokine 
protein levels [20-23] or specific cytokine mRNA levels as the 
readout [24,25]. Human whole blood has also been used as a 
model to assess the whole genome transcriptional response to 
TLR4 [26] or TLR2 and TLR4 ligands [27], although these 
studies only looked at one time point with limited analysis. 

The objective of our study was to undertake a detailed 
comparative analysis of the in vitro global temporal transcriptional 
response to TLR4 and TLR2 ligation in human whole blood. To 
better understand this gene transcriptional response we used a 
variety of bioinformatics approaches to delineate both the 
temporal response and the key transcriptional differences resulting 
from TLR2 and TLR4 ligation and demonstrate that in a whole 
blood system that the response to TLR stimulation can resemble 
that previously identified in isolated immune cells. 

Materials and Methods 

Ethics Statement 

This study was approved by the Central London 3 Research 
Ethics Committee (09/H07 16/41). All participants gave written 
informed consent. 

Human Volunteers 

Six healthy volunteers (self-reported questionnaire); three male, 
three female; aged 25-50 years old; of similar ethnic background 
were recruited into the study. Sixty ml of whole blood from each 
volunteer was taken at 9 am into 10 ml Vacutainers with sodium 
heparin 1 7 international units/ml (ED Vacutainer). 

Whole Blood Cellular Composition 

Measured by CeUtac Automated Hematology Analyzer (MEK- 
6400J/K, Nihon Kohden) at 0 hour, volunteer's results listed 
Table SI. 
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In vitro Whole Blood Stimulation 

In 24 well plates (Costar 3524, Coming Incorporated), 1 ml of 
heparinised whole blood was stimulated either in the presence of a 
final concentration of 200 ng/ml of Pam3CSK4 (Invivogen), 
1 ng/ml of LPS (from Salmonella Minnesota R595, Enzo Life 
Sciences) added in a volume of 100 ^ll with RPMI-1640 with 
GlutaMAX (Life Technologies). Media control samples were 
cultured with the addition of 100 1X1 of RPMI-1640 with 
GlutaMAX. Samples were incubated at ?)T'C., 5% COg for 0, 1. 
3, 6, 12 or 24 hours at which point the contents of the well were 
thoroughly mixed with 2 mis Tempus Solution (Applied Biosys- 
tems/Ambion) to lyse the cells and stabilise the RNA. Samples 
stored at — 80°C until RNA processing. 

Endotoxin Testing 

Reagents (excepting LPS and Pam3CSK4) including Sodium 
heparin Vacutainers were tested for endotoxin contamination by 
Limulus assay and were found to be endotoxin free (<0.03 EU/ 
ml). Pam3CSK4 was tested by the manufacturer and confirmed to 
be endotoxin free (<0.001 EU/p,g). 

RNA Processing 

RNA was isolated using the PerfectPure RNA Blood kit (5- 
PRIME) according to manufacturer's instructions. 2.5 (J.g of 
isolated RNA was globin RNA reduced using the GLOBINclear 
96-weIl format kit (applied Biosystems/Ambion) according to the 
manufacturer's instructions. Isolated and globin reduced RNA 
quantity was assessed using either Nanodrop 1000 or Nanodrop 
8000 spectrophotometer (NanoDrop Products, Thermo Fisher 
Scientific), RNA rjuahty was assessed using an Agilent 2100 
Bioanalyser (Range 6.5-9.5) (Agilent technologies). 200 ng of 
globin reduced RNA was amplified to generate biotinylated 
amplified antisense cRNA using the lUumina CustomPrep RNA 
amplification kit (Applied Biosystems/Ambion). 750 ng of cRNA 
was hybridized to lUumina Human HT-12 V4 BeadChip arrays 
(lUumina) and scanned on lUumina iScan. GenomeStudio 
(Illumina) was used to perform quEility control and generate signal 
intensity. Two samples were excluded from further analysis at this 
stage as they failed quality control measures (0 hour media control 
xl, 6 hour media control xl). 

Microarray Analysis 

Raw background subtracted data was processed using Gene- 
spring V12.6 (Agilent Technologies) and the following principles 
were applied to all analyses. After background subtraction low 
signal values (<10) were then set to a threshold of 10, log2 
transformed and per chip normahsed using 75* percentile shift 
algorithm. Per-transcript normalisation was undertaken by nor- 
malisation to the median of a defined control group. Transcripts 
were then filtered out if they were not significantly (/)<0.01) 
different in intensity value compared to the background in at least 
10% of all the samples. 

The resulting transcripts were then subjected to statistical 
filtering (either One-way ANOVA or 2-way ANOVA) with 
multiple testing correction (Benjamini-Hochberg j&<0.01), fol- 
lowed by a further filtering of transcripts by fold change (FC) in 
which transcripts were filtered if less than 1.8 FC different 
between variables of interest. Expression heatmaps were generated 
within Genespring VI 2.6. Heatmap clustering was undertaken 
using DifiTerential distance metric and Wards linkage rule, unless 
otherwise stated. 

Media controls from 2 volunteers had evidence of activation of 
inflammatory genes by 3 hours of culture (Figure SI A). This 



activation persisted in all of the subsequent time points for these 
individuals (not shown) appears to be revealed in culture, is 
independent of the individual, and independent of the length of 
time in transport conditions (Figure SIB). These samples were 
excluded from the study. 

The data discussed in this publication have been deposited in 
NCBI's Gene Expression Omnibus and are accessible through 
GEO Series accession number GSE55375. 

/c-means Clustering 

Within Genespring 12.6 the normalised significant transcript 
hsts were separately clustered by A-means clustering into 9 clusters 
using Euclidian distance metric. Number of clusters was chosen by 
the number of order branches in the dendrogram from the LPS 
expression heatmap. Clusters were compared across stimulations 
using Pearsons Correlation (within Graphpad Prism V6). 

Ingenuity Pathway Analysis (IPA) 

The canonical pathway, gene function annotation and upstream 
analyses were generated through the usi; of IPA (Ingenuity 
Systems, www.ingenuity.com). Significant transcripts id(;ntified 
from GX microarray analysis were uploaded into IPA. For time 
point analyses these lists were filtered by mean FC (>1.8) 
compared to media control at the time point. 

Type I IFN Regulated List 

A list of human type I regulated genes was obtained from the 
Interferome V2 database (Accessed June 2013) [28]. 

Quantitative PGR 

From the globin reduced RNA cDNA was synthesised using 
High Capacity cDNA Reverse Transcription kit (Applied Biosys- 
tems), according to the manufacturer's instructions followed by 

RNase H (Promega) treatment for 30 min at 37°C. IFNBl, ILIA, 
IL6, NFKBl, NFKB2, STATl, STAT2 and IRF7 gene expression 
were quantified by real-time PCR (79()0HT, Applied Biosystems) 
using the TaqMan system, and normalised to GAPDH mRNA. 

Primer probes used were IFNBl (Hs01077958_sl); ILIA 
(Hs00174092_ml); IL6 (Hs00985639_ml); NFKBl 
(Hs00765730_ml); NFKB2 (Hs010208901); STATl 
(Hs01013996_ml); STAT2 (H.s01013123_ml); IRF7 
(Hs01014809_gl); GAPDH (Hs02758991_gl) (all Applied Biosys- 
tems). 

Results 

Media Controls 

Analysis of the media controls alone over time revealed 377 
significandy expressed transcripts over time. The peak of this 
difference was at 24 hours with 281 transcripts more than 1.8FC 
different compared to the 0 hour samples (Figure S2A). These 
differentially expressed transcripts at 24 hours were enriched for 
inflammatory and metabolic function pathways (Figure S2B). For 
this reasons all data for stimulations are filtered against the media 
controls from the time point of interest and fold changes calculated 
as relative to the media controls at the same time point. 

Identification of Significantly Differentially Regulated 
Transcripts following LPS and Pam3CSK4 Stimulation 

To gain insight into the differential temporal gene expression in 
response to TLR4 and TLR2 ligation we performed a compar- 
ative microarray analysis of LPS and Pam3CSK4 stimulated 
human whole blood and accompanying media controls over a 
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time course. From four healthy human volunteers 1 ml of 
heparinised whole blood was stimulated with either the TLR4 
ligand LPS (1 ng/ml), the TLR2 ligand Pam3CSK4 (200 ng/ml), 
(concentrations shown to result in a plateau in previous studies 
within our laboratory), or incubated only with media as a control. 
RNA was isolated at 0, 1, 3, 6, 12 and 24 hours for each stimulus 
or control. 

LPS (TLR4) stimulation induced a greater number of differen- 
tially regulated transcripts as compared to Pam3CSK4 (TLR2) 
and had a higher magnitude of response. For this reason the two 
stimulations were first analys(;d independently to generate the 
significant transcript lists compared to media control over the time 
course after stimulation. LPS stimulation resulted in the differen- 
tial expression of 4777 transcripts (mapping to 3571 unique genes 
in IPA) whereas Pam3CSK4 stimulation resulted in only 1202 
differentially regulated transcripts (mapping to 922 unique genes 
in IPA) (Figure lA). Expression of these transcripts varied over 
time. Overall 90% of the 1202 significant Pam3CSK4 transcripts 
were shared with the LPS stimulation (Figure IB). These shared 
1093 transcripts when ancdysed by canonical pathway analysis 
(within IPA) were shown to contain TREM, TNF and NFkB 
signalling amongst the top 5 pathways (ranked by significance). 
Canonical Pathway analysis of the 3684 transcripts significandy 
expressed following LPS and not Pam3CSK4 stimulation revealed 
"IFN signalling" as the most significant pathway. 

/c-means Clustering Analysis Reveals Similarities and 
Differences between Stimulations 

A-means clustering was applied to cluster the significant 
transcript lists based on their similarity in expression over time 
for each stimulation separately (fuU composition of the clusters 
listed in EXCEL Files SI & S2). For each cluster the most 

significant canonical pathway (IPA) was determined, reflecting 
gene enrichment within each cluster. The clusters were then 
compared by their expression profile over time (Pearson correla- 
tion, Table S2) and the top canonical pathway to discover similar 
clusters between LPS (termed L) and Pam3CSK4 (termed P). 
Based on these criteria ten clusters in response to both stimulations 
were equivalent. Of the remaining clusters, six shared similar 
expression profiles but had different top canonical pathways and 
two clusters (P4 and P7) had no equivalent in the LPS stimulated 
clusters (Figure 2). 

The similar clusters with "Granulocyte Adhesion and Diape- 
desis" as the top canonical pathway (L7 and P2) were small in 
terms of numbers of genes which were highly expressed by 1 hour 
after stimulation with both LPS and Pam3CSK4 and peaked in 
expression between 3 to 6 hours. 80'/o of the genes in the P2 
cluster were found within the L7 cluster. These common genes 
were predominantiy chemokines and cytokines: C15orf48, CCL2, 
CCL20, CCL3, CCL3L1, CCL3L3, CCL4L1, CCL4L2, CCL8, 
CCRL2, CXCL2, EBB, ILIA, IL6 and TNF. 

The "Acute Phase response signalling" clusters (L4 and P8) 
were characterised by genes involved with inflammatory response. 
Approximately 60% of the genes in the PS cluster were also found 
within the L4 cluster and these common genes included those 
involved with the inflammatory' response: ORMl, ORM2, 
HAMP, IRAK2, PI3, PTGES, TNFAIP6, TNIP3. In addition, 
within the L4 cluster but not the P8 cluster were genes involved 
with IFN regulation: IFI44L, IFITl, IFIT3, IFNG, OAS3, OASL. 

The "Protein Ubiquitination pathway" clusters (LI and PI) had 
104 common genes between stimulations which included heat 
shock protein genes (HSP90AA1, HSP90AB1, HSPA,5, HSPDl, 
HSPEl and HSPHl) and Proteasome PA700/20S genes (PSMAl, 
PSMC3, PSMC4, PSMDl and PSMD14). 



The "Complement system" clusters (in both L8 and P5) were 
characterised by genes that were upregulated, had their highest 
expression at 24 hours and included the metaUothionein genes 
(MTIG, MTIH, MTIE, MTIX, MTIM and MTIF) which were 
amongst the highest expressed genes at 24 hours following both 
LPS and Pam3CSK4 stimulations (Figure S3). Clusters L.t and P6 
contained transcripts that were down regulated over time and 
were similarly enriched for the pathway "Dermatan Sulfate 
Biosynthesis (Late stages)". 

Cluster L2, which peaked at 6 hours and was characterised by 
"Activation of IRF by cytosolic pattern recognition rec:eptors" as 
the top pathway did not ha\c a similar pathway in the Pam3CSK4 
clusters. Examination of tiu- genes within this pathway revealed 
that there were detected a number of IFN regulated/regulatory 
genes (DDX58, DHX58, IFIHl, IFIT2, IFNBl, IRF7, STAT2). 

Clusters L3, LO, L6 and P4, P7 and PO were all enriched for 
immune function pathways although there were no similarities 
between the stimulations in terms of top significant pathway and 
kinetic profile. Cluster P3 was enriched for the pathway "UDP-N- 
acetyl-D-galactosamine Biosythesis 11" and there was no simflar 
cluster following LPS stimulation. 

Common and Distinct Changes in Transcription at the 
Different Time Points 

To better understand the temporal response we undertook 
analysis at each time point using the significant transcript lists. 
Significantly expressed transcripts at each time point were 
identified by comparing the mean expression of the differentially 
regulated transcripts in response to LPS and Pam3CSK4 
stimulation at each time point to the media control at that time 
point and filtering those transcripts which were less than 1.8 FC 
different to the media control (full listings of identified transcripts 
given as EXCEL Files S3 & S4). From this it was observed that the 
peak of the transcriptional response compared to media control 
occurred at 6 hours for both LPS and Pam3CSK4 stimulations 
(Figure 3A). This diflerential transcriptional response was still 
evident at 24 hours foUowing in vitro LPS stimulation of whole 
blood. 

At each time point the per-time point transcript lists were 
analysed by canonical pathway analysis within IPA. The top 25 
significant pathways (by mean -log p value of the LPS time point 
pathways and compared to Pam3CSK4) are shown (Figure 3B). 
This analysis revealed a large number of canonical pathways with 
simflar kinetics of significance in both LPS and Pam3CSK4 
stimulations, with the exception of the "IFN Signalling" pathway 
where there was a clear difierence in significance between 
stimulations. 

The pathways that were most significant at 1 hour and then 
diminished in significance over time were "Agranulocyte Adhesion 
and Diapedesis", "Granulocyte Adhesion and Diapedesis", 
"Diflerential Regulation of cytokines production in Epithelial 
Cells by IL-17A and IL-17F", "Differential regulation of cytokines 
production in macrophages and T-Helper CeUs by IL-17A" and 
the pathway "Role of Hypercytokinemia/hyperchemkinemia". AU 
the genes contributing to these pathways at this time point were 
chemokine and cytokine genes. "IL-10 signalling", "Communica- 
tion between Innate and Adaptive Immune ceUs" and "Athero- 
sclerosis Signalling" pathways were also most significant at 1 hour 
and then diminished in significance over time. In addition to 
cytokine and chemokines genes the "Atherosclerosis SignaUing" 
and "Communication between Innate and Adaptive Immune 
cells" pathways had CD40, ICAMl, ORMl and ORM2 as being 
significandy expressed within the pathways at 1 hour following 
LPS and Pam3CSK4 stimulation. The "IL-10 signalling" pathway 
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Figure 1. LPS or Pam3CSK4 stimulations results in a differential response in gene expression overtime. (A) 1 ml of human whole blood 
from healthy volunteers (N = 4) was stimulated with either Pam3CSK4 (200 ng/ml), LPS (1 ng/ml) or media control for different lengths of time (0, 1, 3, 
6, 12 and 24 hours). Stimulations were analysed independently: media control compared to Pam3CSK4 and media control compared to LPS revealed 
1202 and 4777 significantly expressed transcripts respectively. Transcripts were identified by normalising expression values to the median of the 0 
hour samples, filtering by detection from background, statistical filtering (2 way ANOVA with Benjamini Hochberg multiple testing correction p<0.01 ) 
and retaining transcripts whose expression was greater than 1.8 PC different between the media control and stimulation samples at one or more 
time point. (B) A Venn diagram of both significant transcript lists. Within the Venn for each subset the number of transcripts is given, with unique 
genes within IPA in brackets. For transcript lists the top 5 canonical pathways (IPA) are shown as well as a heat map of the normalised expression 
values of these transcripts for both stimulations over time. 
doi:1 0.1 371/journal.pone.0097702.g001 



in addition to chemokiiies and cytokines liad tlie transcription 
factors, JUN, NFKB I, NFKB2, NFKBIA, NFKBIE as well as tlie 
gene SOCS3 within the patliway. 

The pathways which were most significant at 3 hours following 
both LPS and Pam3CSK4 stimulation were "NFkB Signalling" 
and "IL-6 Signalling". There were a number of pathways that 
were significant by 1 hour and remained significant following both 
LPS and Pam3CSK4 stimulations: "TNFR2 signalling", 
"TREMl signaUing", Altered T ceU and B Cell signaUing in 
Rheumatoid arthritis" and "Dendritic cell maturation" pathways. 

The "IFN signalling" pathway in response to LPS was only 
significant (/)<0.01) from 3 hours onwards and was most 
significant at 6 hours (Figure 3B). There was a clear difference 
between the level of significance of the IFN signalling pathway 



following LPS and Pam3CSK4 stimulation, with the "IFN 
signalling" pathway not reaching the threshold of significance 
(/)<0.01) at any time point following Pam3CSK4 stimulation. 

The significantly expressed genes at 1 hour were analysed, as 
from the canonical pathway analysis we had identified enrichment 
of cytokines and chemokines genes at this time point. In addition 
this was the earliest measured time point following stimulation, 
when autocrine and paracrine signalling leading to induction of 
mRNA should be at its minimum and, 334 LPS and 165 
Pam3CSK4 genes respectively were more than 1.8 FC different to 
media control at 1 hour. Cytokines/ chemokines and transcrip- 
tional regxilators (identified by IPA gene function annotation) 
combined accounted for approximately 20% of the significantly 
expressed genes at this time point. There was a large degree of 
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Figure 2. Ar-means clustering of the significant transcript lists reveals similar clusters. Mean normalised expression profile of individual 
clusters (y axis, ±SD, n = 4) at each time point (x axis; 0, 1, 3, 6, 12 and 24 hours) for each cluster. N: number of transcripts within cluster; P: most 
significant canonical pathway (IPA). Clusters are grouped by similarity in kinetic profile (Pearsons correlation. Table S2) and top canonical pathway. 
doi:1 0.1 371/journal.pone.0097702.g002 



overlap between stimulations of the significant cytokines/ chemo- 
kines at 1 hour, with all of the 17 significant cytokines/chemokines 
genes from the Pam3CSK4 list also being identified in the LPS list 
(Figure 3C). Following LPS but not Pam3CSK4 stimulation, of 
note IL12B, IFNBl, IFNG, IL23A and CXCLIO were also 
significandy expressed as early as 1 hour (Figure 3C). There were a 
number of significandy expressed transcriptional regulators by 1 
hour, with the majority of significandy expressed transcriptional 
regulators identified in the Pam3CSK4 1 hour list shared with 
LPS, including tiiose of the NFkB famUy (NFKBl, NFKB2, 
NFKBIA, NFKBIE, NFKBIZ and REL) and the AP-l/CREB 
family JUN, JUNE, ATF3, BATF3). Transcriptional regulators 
involved in cell development, proliferation and differentiation 
(BTG2, EGR3, ETV3, HHEX, MAFF, SKIL) and the post- 
transcriptional regulator ZFP36 (also known as TristetraproUn) 
were also identified (Figure 3C). 

Upstream Analysis Identifies Potential Transcriptional 
Regulators 

In order to identify which transcriptional regulators may be 
responsible for the differences obser\'ed in gene expression 
between LPS and Pam3CSK4 stimulation we undertook upstream 
analysis of the gene lists at each time point within IPA. Upstream 
analysis attempts to predict which transcriptional regulator may be 
responsible for the observed differential gene expression by 
comparing the genes known to be regulated by a transcriptional 
regulator (derived from the literature) to those significandy 
differentially expressed genes identified at each time point from 
this analysis. It can he seen that the top predicted transcription 
regulator for both stimulations was the NFkB complex and that it 
was predicted to be activated from 1 hour onwards (Figure 4A). 

By comparison of the mean mRNA expression of the predicted 
transcriptional regulators compared to media control it can be 
seen that the magnitude and temporal response was similar for 
LPS and Pam3CSK4 stimulation for PPRCl and SPl (Figure 4B). 
The predicted transcriptional regulators NFkBI, CEBPA, EGRl 
shared similar mRNA temporal profiles between stimulations; 
however the magnitude of the response was greater upon LPS 
stimulation. RELA, EZH2 and STAT3 were predicted to be 
activated in both stimulations however had different temporal 
profiles between stimulations. RELA, EZH2 and STAT3 were 
upregulated earlier and to a greater magnitude in LPS compared 
to Pam3CSK4 and this upregulation persisted over time in the 
LPS stimulation compared to Pam3CSK4 stimulation. The 
mRNA expression of STATl and IRFl, which were only 
predicted to be activated following LPS stimulation, had different 
profiles being induced earlier and to a greater magnitude in LPS 
compared to Pam3CSK4 stimulations. 

The NFkB genes REL, RELA, RELB, NFkBI and NFkB2 
were all significandy expressed and present in both LPS 4777 
transcript and Pam3CSK4 1202 transcript lists (Figure 40). 
Transcriptional r(;gulators imoK'cd in IFN regulation IRFl, IRF3 
and STATl were predicted to be ac:tivated in the LPS stimulation 
but not the Pam3CSK4 stimulation (Figure 4A). IRFl, IRF2, 
IRF4, IRF7, IRF8, IRF9, STATl, STAT2, STAT3, STAT4 and 
STATS A were aU identified in the 4777 LPS and not the 
Pam3CSK4 1202 transcript list as being significantly differentially 
expressed, IRF7 having the greatest FC induction compared to 
media control. IRF3, IRF5 and IRF6 expression was not 



significandy regulated following LPS or Pam3CSK4 stimulation 
compared to media control (Figure 4D). Overall the inducible IRF 
and STAT gene expression level peaked later (6 hours) in 
expression at 6 hours compared to that of the NFkB family of 
genes which peaked between 3 to 6 hours (Figure 4E). We 
validated this difference between LPS and Pam3CSK4 in terms of 
the magnitude of expression at the 6 hour time point for NFKBI, 
NFKB2, STATl, STAT2 and IRF7 by quantitative real time 
PGR (Figure S5). In addition the FG relative to the media control 
was lower for STATl, STAT2 and IRF7 compared to tiie NFKBI 
and NFKB2 foUowing Pam3CSK4 stimulation. 

Upstream Analysis Identifies Potentially Active Cytokines 

In order to attempt to identify the cytokines that were 
potentially involved in autocrine /paracrine signalling and subse- 
quent gene cxprc-ssion regulation in response to LPS or 
Pam3CSK4 stimulation, we again undertook per-time point 
upstream analysis within IPA. TNF, ILIB and ILIA were 
predicted as potentially activated cytokines, and this predicted 
activation was early (by 1 hour) and sustained (Figure 5A). The 
mRNA expression of these cytokines compared to media control 
revealed them to be highly upregulated by 1 hour in both LPS and 
Pam3CSK4 stimulations. Although the magnitude was higher in 
LPS stimulation, they shared similar kinetic profile between 
stimulations (Figure 5B). 

IL12B was predicted to be activated upon LPS stimulation but 
not Pam3CSK4 and this difference was observed in the mean 
mRNA expression of the IL12B gene (encoding for IL12p40) 
which was upregulated in response to LPS compared to media 
control and not following Pam3CSK4 stimulation (Figure 5). IFN 
cytokine activation was also only predicted in the LPS stimulation 
(IFNA2, IFNLl and IFNBl), in keeping with this the kinetic 
profile of IFNBl can be seen to be different between LPS and 
Pam3CSK4 stimulations. This difference between LPS and 
Pam3CSK4 stimulation was validated by real time PGR at the 
peak of IFNBl expression at 3 hours. ILIA and IL6 which also 
peaked at 3 hours were seen to have a robust response following 
both LPS and Pam3GSK4 stimulation, although with a much 
greater response following LPS (Figure S5). IFNA and IFNL genes 
were of low magnitude expression (< 1 .8 FG different to media 
control) under both conditions of stimulations (not shown). 

IFN Gene Expression is Dominant Following LPS but not 
Pam3CSK4 Stimulation 

Our data show that as early as 1 hour post stimulation IFN gene 
expression is seen to be upregulated foUowing LPS but not 
Pam3GSK4 stimulation. This difference in IFN signalling was 
further emphasised by the A:-means clustering, canonical pathway 
analysis and upstream analysis of both potential transcriptional 
regulators and cytokines highlighting a difference in IFN signalling 
following LPS and Pam3CSK4 stimulation. 

IRF and STAT genes had been identified by k-me.ans clustering, 
predicted upstream analysis, and canonical pathway analysis as 
being differentially activated between LPS and Pam3CSK4 
stimulations. 

To test if this difference in expression of IRF and STAT genes 
resulted in differential expression of IFN regulated genes we used a 
list of human type 1 IFN regulated genes generated from the 
Interferome database v2.0. We compared the expression of these 
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Figure 3. Transcript lists analysed at each time point. (A) A graph showing the number of genes from the respective significant transcript lists 
(4777 LPS and 1202 Pam3CSK4 lists) at each time point which are more than 1.8 FC different compared to the media control at that time point. (B) 
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The significantly expressed transcript lists (1202 Pam3CSK4 transcript list and 4777 LPS transcript list) were analysed in IPA. For each time point only 
genes whose expression were 1.8 FC different from the media control at that time point were taken into consideration. Shown is a heatmap of 
pathway significance of the top 25 IPA canonical pathways for each time point where significance criteria met (Fishers Exact test p<0.01). The IPA 
canonical pathways were chosen by identifying from the LPS stimulation analyses the top 25 most significant pathways across the time points (mean 
-log p value) and then compared to Pam3CSK4. (C) Venn diagrams of cytokine/chemokines and transcriptional regulators identified using IPA gene 
functional classification from LPS 4777 and Pam3CSK4 1202 transcript lists with mean expression greater than 1.8 FC different to media control at 1 
hour. Listed adjacent to the Venn diagrams are the genes from each subset. 
doi:1 0.1 371/journal.pone.0097702.g003 



type 1 IFN regulated genes between LPS and Pam3CSK4 
stimulations. From this it can be seen that LPS stimulation results 
in a greater number of difTerentially regulated type 1 IFN genes 
compared to Pam3CSK4. In addition the magnitude of this 
differential regulation was much higher in LPS stimulated whole 
blood compared to Pam3CSK4 stimulated whole blood and the 
peak of this transcriptional response following both LPS and 
Pam3CSK4 was at 6 hours. (Figure S4). 

Induction of NFkB Signalling Pathway is Similar 
Following LPS and Pam3CSK4 Stimulation in Contrast to 
the IFN Signalling Pathway 

The similarity in NFkB signalling and difference in IFN 
signalling following Pam3CSK4 and LPS stimulation is reflected 
in the IPA canonical pathways. There is similarit}' in the IPA 
canonical pathway "NFkB signalling" at 3 hours (the peak of 
significance for this pathway following LPS and Pam3CSK4 
stimulations) following LPS and Pam3CSK4 stimulation. However 
there is clear difference IFN Signalling canonical pathway at 6 
hours (the peak of significance for this pathway for LPS, the 
pathway is not significant at any time point following Pam3CSK4 
stimulation) where there is seen to be difference following LPS and 
Pam3CSK4 stimulations (Figure 6). 

Comparison of Differential Transcriptional Expression 
Following Either in vitro Stimulation or in vivo LPS 
Stimulation of Human Whole Blood 

In vivo LPS stimulated human whole blood data was obtained 
from the Gene Expression Omnibus (GSE3824 [16]), we analysed 
this data within Genespring (vl2.6) using the same methodological 
principles used to analyse our own data. It can be seen that as 
opposed to in vitro LPS stimulation, following human in vivo LPS 
administration the transcriptional response had returned to 
baseline by 24 hours (Figure 7A). Comparative analysis of the 
significantly differentially regulated gene lists over 24 hours (per- 
time point analysis was not possible due to the different time points 
sampled, transcript level analysis not appropriate due to the 
difference in technology platforms), rcveak-d that approximately 
40% of the significantiy differentially regulated genes following 
in vivo LPS administration were also differentially regulated 
following i?i vitro LPS stimulation. Canonical pathway analysis 
revealed these shared genes between stimulations were enriched 
for IFN, TNFRl, TREM and NFkB Signalling. 

The genes differentially regulated following in vivo and not 
in vitro LPS stimulation were enriched for pathways involved in cell 
cycle, cellular development as well as immune function, however 
the most significant pathways for those genes differentially 
regulated following in vitro and not in vivo LPS stimulation were 
predominantly immune function pathways. 

Transcriptional regulators which were differentially expressed 
following both in vivo and in vitro LPS stimulation included NFkB 
family genes (RELA, RELB, NFKBl, NKFB2) and IRF and 
STAT transcriptional regulators (IRFl, IRF2, IRF7, IRF8, IRF9, 
STATl, STAT2, STAT3, STAT4, STAT5A) (Figure 7B). Anal- 
ysis of the kinetic profile of these transcriptional regulators showed 



similar timing of the peak of induction for both in vivo and in vitro. 
Of note in-vivo expression then diminished soon after the peak of 
expression in contrast to that observed in vitro where persistent 
differential expression was observed stiU at 24 hours. 

ILIA, ILIB, IL15, IL32 and TNF were cytokines which had 
been identified by in vitro upstream analysis and seen to be 
differentially expressed in-vitm and were also significantly differen- 
tially expressed following in vivo LPS administration (Figure 7C). 
However for ILIA, ILIB and TNF it can be seen that the 
magnitude of this induction is much lower in vivo than in vitro and 
that the kinetic profiles are different with expression returning to 
that of the media controls by 9 hours following stimulation in vivo 
in contrast to the in vitro (^xpc-riment where expression remains 
elevated for the duration of the experiment compared to media 
controls. 

The following cytokines/ chemokines which had been predicted 
to be activated cytokines from the in vitro experiment and which 
were significantly differentially expressed in vitro but not in vivo 
were IFNBl, IFNG, CSF2, IL5, IL6, ILIO, IL12B, and IL21. 

Discussion 

A central goal of our study was to identify the key transcrip- 
tional differences between TLR4 and TLR2 ligation in a human 
whole blood system and demonstrate that in a whole blood system 
that the response to TLR stimulation can resemble that previously 
identified in isolated immune cells. We show that the magnitude of 
the transcriptional response both in terms of the number of 
differentially regulated genes as well as level of mRNA expression 
following LPS stimulation was greater and more sustained than 
the transcriptional response to Pam3CSK4. There was a common 
transcriptional response following stimulation with LPS and 
Pam3CSK4 which was dominated by NFkB regulated genes. In 
addition there was a separate IFN regulated transcriptional 
response (IRF and STAT mediated) seen mainly following LPS 
stimulation. This significant difference in the IFN response 
between TLR4 and TLR2 ligation could be seen as early as 1 
hour post stimulation in our study, and reinforced in the later time 
points despite the potentially complex autocrine and paracrine 
signalling in the human whole blood system. 

The common early transcriptional response following both LPS 
and Pam3CSK4 stimulations at 1 hour was characterised by 
highly expressed cytokines and chemokines (CCL2, CCL20, 
CCL23, CCL3, CCL3L1/CCCL3L3, CCL4, CXCLl, CXCL2, 
EBI3, EDNl, ILIA, ILIB, ILIRN, IL6, IL8 and TNF) as well as 
the significant upregulation of transcriptional regulators including 
tiie NFkB family (NFKBl, NFKB2, NFKBIA, NFKBIE, NFKBIZ 
and REL) and AP-l/CREB (JUN, JUNB, ATF3, BATF3). The 
gene ZFP36 which encodes for tlu; j)r()t(;in Tristetraprolin was also 
identified as significantiy up regulated by 1 hour which has been 
shown to act in a post-transcriptional rc-gulatory role by binding to 
the mRNA of some cytokines and promoting their degradation 
[29]. The peak of expression compared to media control of NFkB 
genes was 3 to 6 hours following both LPS and Pam3CSK4 
stimulation. By 24 hours following Pam3CSK4 stimulation the 
expression of these NFkB genes was not significantiy different 
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Figure 4. Predicted transcriptional regulator identification and NFkB, IRF and STAT gene expression following LPS and Pam3CSK4 
stimulation. (A) Predicted upstream transcriptional regulators from IPA; stimulations analysed independently using Pam3CSK4 1202 and LPS 4777 
transcripts lists, at each time point only genes whose expression were 1.8 PC different from the media control at that time point were taken into 
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consideration. Predicted upstream transcription regulators which met the criteria (p<0.01 (Fishers Exact Test) and z activation score >2.5) shown 
plotted by z-activation score only at the time points where significance criteria met. (B) Mean mRNA expression of predicted transcription regulators 
plotted as log2 fold change (y axis) across time (x axis; 0, 1, 3, 6, 12 and 24 hours), fold change is relative to the media control at each time point, only 
those predicted transcription regulators whose mRNA expression is >1.8 FC relative to media control at one or more one time points are shown. (C) 
NFkB genes. IVIean mRNA expression of the NFkB family genes. All 5 genes were significantly expressed and present in both Pam3CSK4 1202 and 
LPS 4777 transcript lists. (D) Interferon Regulatory Factors. Mean mRNA expression of the IRF and STAT genes. IRF1, IRF4, IRF7, IRF8, IRF9, ST ATI, 
STAT2, STAT3, STAT4 and STAT5A were present in the LPS 4777 significantly expressed transcript list, none of the IRF or STAT genes (except STAT5B) 
were present in the 1202 Pam3CSK4 significant transcript list. (E) Temporal Kinetics of NFkB and induced IRF and STAT genes. Plotted for 
both LPS and Pam3CSK4are mean fold change relative to media controls of the NFkB genes (NFKB1, NFKB2, REL, RELA, RELB) and mean fold change 
relative to media controls of selected IRF genes (IRF1, IRF2, IRF4, IRF7, IRF8, IRF9, STAT1, STAT2, STATS, STAT4 and STAT5A). 
doi:10.1371/journal.pone.0097702.g004 



compared to media control, and similarly the expression of many 
proinflammatory cytokines (e.g. ILIA, IFNG and TNF) had 
diminished significantly by 24 hours compared to media control. 
This is in contrast to the expression following LPS stimulation 
where these NFkB genes and cytokines were still significandy 



upregulated compared to media control at 24 hours. This 
prolongation of the NFkB signalling seen following LPS stimula- 
tion maybe as a result of secondary induction of NFkB via the 
TRIF-TRAM adaptor molecules, or as a consequence of the IFN 
signalling [7,8] . We also identified a common group of metaUo- 
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Figure 5. Identifying potential cytokines involved in autocrine gene regulation by upstream analysis within IPA. (A) Predicted 
activated cytokines from IPA upstream analysis; stimulations analysed independently using Pam3CSK4 1202 and LPS 4777 transcripts lists, at each 
time point only genes whose expression were 1.8 fold different from the media control at that time point were taken into consideration. Predicted 
upstream cytokines which met the criteria (p value <0.01 (Fishers Exact Test) and z- activation score >2.5) shown plotted by z-activation score only at 
the time points where significance criteria met. (B) Cytokine identified from either predicted upstream analysis or canonical pathway analysis (Fig. 3) 
mean mRNA expression plotted as log2 fold change (y axis) across time (x axis; 0, 1, 3, 6, 12 and 24 hours), fold change is relative to the media control 
at each time point, only those predicted cytokines whose mRNA expression is >1.8 fold upregulated relative to media control at one or more one 
time points are shown. 
doi:1 0.1 371/journal.pone.0097702.g005 
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Figure 6. NFkB and Interferon signalling pathways. Shown at the peak of their significance, 3 and 6 hours respectively. Significantly expressed 
genes within the pathway (from Pam3CSK4 1202 and LPS 4777 lists) shaded red if upregulated or blue if down regulated. 
doi:1 0.1 371/journal.pone.0097702.g006 



thionein genes (MTIG, MTIH, MTIE, MTIX, MTIM, MTIF) 
as being amongst the most up regulated genes at 24 hours 
following both LPS and Pam3CSK4 stimulation. Metallothionein 
genes are known to have a key role in oxidative stress and heavy 
metal detoxification [30] . It has also been shown that they may be 
induced by cytokines such as IL-6 and they are thought to have a 
gene regulatory role in inflammation [31-33]. The late induction 
of these genes in our system suggests that they may be being 
induced by subsequent auto/paracrine signalling downstream of 
the MyD88 adaptor. 

Following LPS stimulation a differential transcriptional response 
of IFN regulated/ regulatory genes was observed as compared to 
Pam3CSK4. This difference was detectable as early as 1 hour post 
LPS stimulation with the significant upregulation of the cytokine 
genes IFNBl and IFNG and the transcriptional regulator IRF8 
following LPS and not Pam3CSK4 stimulation. The IRF and 
STAT (IRFl, IRF2, IRF4, IRF7, IRF8, IRF9, STATl, STAT2, 
STAT3, STAT4 and STAT5A) genes were seen to be differen- 
tially regulated following LPS stimulation and the peak of their 
expression was at 6 hours. This difference in induction of IRF and 
STAT transcription regulators following TLR4 and not TLR2 
ligation is consistent with previous studies which have compared 
the temporal transcriptional response following TLR2 and TLR4 
ligation in murine DCs and macrophages [10,11]. The IFN 



signalling induced following TLR4 ligation and not TLR2 ligation 
can be accounted for by signalling through the MyD88 
independent TRIF-TRAM, IRF3 pathway foUowing TLR4 
ligation resulting in induction of IFNfi, which then leads to 
positive feedback of IFN gene induction through induction of 
IRF7 amongst other IRFs [8] . IRF3 is thought to be constitutively 
expressed and its expression is not thought to be induced by TLR 
ligation, or type I or II IFN. This may explain why IRF3 was 
predicted by upstream analysis to be activated and its mRNA 
expression not significandy induced following TLR4 or TLR2 
ligation in our study [8]. IRF7, which is known to be induced by 
type I IFN subsequent to IRF3 activation by TLR4 ligation was 
seen to be the most induced of the IRFs following TLR4 ligation 
[8] . This stepwise induction of the IRFs may account for the later 
peak of expression at 6 hours of the inducible IRF and STAT 
genes compared to the NFkB genes which peaked at 3 hours in 
our study. This later peak in transcription factors associated with 
IFN regulation compared to NFkB genes is in agreement with 
studies in murine macrophages as well as ex-vivo human blood 
leukocytes. Slight differences in the timing of peaks observed in our 
study and these other studies could potentially be explained by the 
difference in time points sampled or differences in the systems used 
[10,34]. 
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Figure 7. Comparison of WKiole blood in vitro LPS stimulation and in vivo LPS stimulated blood leukocytes. GSE3284 data obtained 
from Gene Expression Omnibus (GEO) repository. Stimulations analysed independently to generate significant transcript lists: 4777 and 4791 
significantly expressed transcripts respectively. Transcripts identified by normalising expression values to the median of the 0 hour samples, filtering 
by detection from background, statistical filtering (2 way ANOVA with Benjamini Hochberg multiple testing correction p<0.01) and retaining 
transcripts whose expression was greater than 1.8 PC different between the media control and stimulation samples at one or more time point. 
Heatmaps of significant transcript lists shown. Venn diagram of unique genes lists with canonical pathway analysis within IPA of subgroups (top ten 
significant pathways shown) (B) Temporal Kinetics of NFkB and induced IRF and STAT genes. For both in vivo and in vitro LPS stimulations 
mean fold change relative to controls of the NFkB genes (NFKB1, NFKB2, REL, RELA, RELB) and mean fold change relative to media controls of selected 
IRF genes (IRF1, IRF2, IRF4, IRF7, IRF8, IRF9, STAT1, STAT2, STAT3, STAT4 and STAT5A). (C) Comparison of cytokine expression. IVlean mRNA 
expression of cytokines which were identified as being significant from the in vivo upstream analysis and were observed to be significantly expressed 
in both in vivo and in vitro datasets were plotted as log2 fold change (y axis) across time, fold change is relative to their controls at each time point. 
doi:1 0.1 371/journal.pone.0097702.g007 



Following LPS stimulation the peak of differential transcription 
compared to media control was at 6 hours in terms of the number 
of differentially regulated genes. At 24 hours 1443 genes were still 
significantly differentially regulated. This persistence of differential 
transcription at 24 hours is observed in human LPS stimulated 
in vitro transcriptional studies [12,13]. H()vve\'er this in contrast to 
that noted in human in-vivo stimulation where transcriptional 
difference peaked between 4 to 6 hours and by 24 hours the 
transcriptional signature had returned to baseline [16]. Several 
cytokines /chemokines were not seen to be significantly induced 
in vivo compared to in vitro and these differences in expression and 
persistence of differential expression observed in vitro compared to 
in vivo could be explained by trafficking/removal of activated 
immune cells out of the circulation in vivo, which is not possible to 
be represented in an in vitro whole blood system. 

The ability to utilise in vitro stimulated human whole blood for 
transcriptomic analysis of the early innate immune response has 
potential advantages over the use of isolated immune cells as whole 
blood stimulation can be carried out in a laboratory where the 
expertise or equipment to isolate immune cells from blood is 
lacking. In addition the volume of blood needed for stimulation is 
much less than that required to isolate immune cells, meaning 
experiments can be performed in populations where access to 
larger volumes of blood is not possible e.g. paediatrics or it could 
be possible to undertake more stimulations/time points with a 
given volume. The results however must be interpreted in the 
context of the complex autocrine/paracrine interactions which 
occur in a mixed cell culture and therefore stimulation of isolated 
immune cells wiU remain advantageous to interrogate in detail a 
specific transcriptional response. 

We have shown that human whole blood can be used to study 
the early temporal transcriptional response following TLR2 and 
TLR4 ligation and that the results obtained are comparable to 
those from isolated murine and human immune cells. 

Supporting Information 

Figure SI Activation of samples by culture conditions. 

(A) Heatmap of expression, clustered by transcripts, shows that 2 
individuals (out of 6) media controls show activation (marked with 
arrows), the genes differentially expressed in these 2 media control 
sampk-s arc similar to LPS, but lower in magnitude. Transcripts 
identified by normalisation to 0 hour samples, filtering by 
detection from background, statistical filtering (2 way ANOVA 
with Benjamini Hochberg p<0.01) and then transcripts retained 
whose expression was >1.8 FC from another condition. (B) 
Activation is not due to length of time in transport conditions, nor 
is it individual specific. The single individual shown here had not 
previously activated, and was included in final dataset. Blood was 
collected at time point 0 and left in sealed vacutainers, the 
vacutainers were opened at one hour intervals and 1 ml human 
whole blood was either immediately mixed with Tempus solution 



(labelled as Direct from vacutainer) or plated for 3 hours with 
either media control (RPMI-1640 with GlutaMAX) or LPS (1 ng/ 
ml) and then mixed with Tempus solution. Reagents and 
containers (including vacutainers) are endotoxin free (undetectable 
by Limulus assay - sensitivity <0.03 EU/ml). Heatmap of 
expression (duplicate stimulations from the same individual shown, 
2619 transcripts), clustered by transcripts shows that regardless of 
length of time in vacutainer activation occurred in all media 
control samples, and is not observed in the direct from vacutainer 
samples, implying that activation is dependent on culture 
conditions and is not a function of length of time ex-vivo or length 
of time spent in the vacutainers. Transcripts identified by 
normalisation to median of "Direi:t frcjm vacutainer" samples, 
filtering by detection from background, statistical filtering 
(ANOVA with Benjamini Hochberg p<Q.Q\) and then transcripts 
retained whose expression was >1.8 FC from another condition. 
(TIF) 

Figure S2 Transcriptional changes in media controls 
over time. (A) Heatmap of normalised expression values of 377 
transcripts which were identified to be significantly differentially 
expressed overtime (transcripts identified by normalisation to 0 
hour samples, filtering by detection from background, statistical 
filtering (One way ANOVA with Benjamini Hochberg p<0.01) 
and then transcripts retained whose expression was > 1 .8 FC from 
the 0 hour samples. Plotted above the heatmaps is the number of 
signiiicantiy expressed genes (mapped in IPA) that were differen- 
tially expressed at each time point by more than 1 .8 FC compared 
to the 0 hour samples. (B) Top ten IPA canonical pathways of the 
significantiy expressed genes at 24 hours, with the -log p \'alue for 
the pathway and the significantiy differentially expressed genes 
listed for each pathway. 
(TIF) 

Figure S3 Metallothionein gene expression. (A) Heatmap 
of averaged Metallothionein mRNA expression over time 
following LPS or Pam3CSK4 stimulation, values normsdised to 
the median of the 0 hour. Note asynchronous scale. 
(TIF) 

Figure S4 Interferon regulated genes. Heatmap of aver- 
aged expression values of Type 1 Interferon regulated genes (List 

obtained from Interferome v2.0), normalised to the median of the 
0 hour, genes retained if they were expressed greater than 1 .8 FC 
from media control in at least one stimulation in one or more time 
points (resulting in 1105 genes). Graphed above heatmap is the 
mean absolute fold change of these Type 1 interferon regulated 
genes. 
(TIF) 

Figure S5 Real time PGR. Real time PCR of selected genes 
following LPS and Pam3CSK4 stimulations and media controls, 
normalised to GAPDH expression. Mean fold change calculated 
between media controls and stimulations. 
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(TIF) 

Table SI Volunteer whole blood composition measured 
by Celltac Automated Hematology Analyzer (MEK- 
64O0J/K, Nihon Kohden) at time point 0 hour. 

(TIF) 

Table S2 Pearson correlations for k-means derived 
clusters from Figure 2. 

(TIF) 

File SI Listings of transcripts, LPS k-means clusters 
from Figure 2. 

(XLSX) 

File S2 Listings of transcripts, Pam3CSK4 k-means 
clusters from Figure 2. 

(XLSX) 

File S3 LPS time course data. Listings of transcripts from 
4777 significant transcript list whose mean expression was 1.8 FC 
different to media control at each time point from Figure 3. 

(XLSX) 



File S4 Pam3CSK4 time course data. Listings of transcripts 
from 1202 significant transcript list whose mean expression was 
1 .8 FC difiFerent to media control at each time point from Figure 3. 
PCLSX) 
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